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We review single-qubit quantum process tomography for trace-preserving and nontrace-preserving 
processes, and derive explicit forms of the general constraints for htting experimental data. These 
forms provide additional insight into the structure of the process matrix as well as reveal a tighter 
bound on the trace of a nontrace-preserving process than has been previously stated. We also 
describe, for completeness, how to incorporate measured imperfect input states. 


Quantum process tomography (QPT) is the gold stan¬ 
dard method to measure and quantify how a quantum 
device-under-test (DUT), such as a quantum gate or a 
transmission channel, will transform some arbitrary in¬ 
put quantum state [HE]. Often a quantum process is 
called a map, as it maps any input state to an out¬ 
put state reflecting the input state’s interaction with 
the DUT. In this note focused on single qubit QPT, 
we review the treatment of a trace-preserving process 
(e.g., a wave-plate) and the treatment of a nontrace¬ 
preserving process (i.e., there are qubit loss or leakage 
errors, e.g., the impact of a polarizer on optical qubits 
may cause qubits to be absorbed), and extract explicit, 
implementable expressions from the general constraints 
given in DEj. These explicit forms, which involve dif¬ 
ferent elements of the process matrix, then provide in¬ 
sight into the structure of the process matrix, and thus 
can serve as useful tools for an experimentalist interested 
in measuring quantum gates with error models including 
qubit leakage. A tighter bound on the trace of a nontrace¬ 
preserving process also emerges, which we point out. 

Furthermore, we discuss input quantum states that are 
not perfectly prepared, which is important to consider 
to reduce experimental errors in the determined process 
matrix. Consequently, we provide a slight modification 
of the methodology for single qubit tomography |5] to 
include such imperfectly prepared input states. To carry 
out single qubit QPT experimentally, one must prepare 
four single-qubit states which will be used to probe the 
DUT. Even if states are prepared with high fidelity, they 
are not perfect, and as such depart from the assump¬ 
tion made in [2]. In practice, one may characterize the 
probe states via quantum state tomography (QST) [31- 
E]. By characterizing the input probe states with QST, 
the quantum process may be more accurately assessed. 
While QST can itself introduce errors, a well calibrated 
QST system minimizes such errors. 
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Device-Under-Test (DUT) 


FIG. 1. The input state p changes to p' upon traversal 
through a device-under-test (DUT). 


After the probes emerge from the DUT, they are mea¬ 
sured, resulting in a total of eight quantum states (four 
input and four output), which are measured by QST. 
From the measured density matrices, the mapping of the 
DUT is computed. The map is normally described as a 
4x4 process matrix operator (y). 

Fig.[T] shows a DUT upon which qubits impinge cor¬ 
responding to a quantum state described by the density 
matrix p. The output qubits’ density matrix is denoted 
by p'. Ordinarily QST produces normailized states; how¬ 
ever the measurements contain additional information on 
the loss. To use the loss information, the density matrix 
of the output state p’ includes a scaling factor (< 1) to 
account for any loss of qubits as they traverse the DUT, 
that is to say, they are not normalized to the output 
qubit flux but to the input qubit flux. In other words, 
while Tr{p) = 1 always, Tr{p') < 1, with the inequality 
holding for a non-trace-preserving process. For example. 


P' = 


1/2 O' 
0 0 


1/2 


1 0 
0 0 


( 1 ) 


has trace equal to | and may be interpreted as an ensem¬ 
ble of qubits in the pure state |0), reduced in quantity by 
50%. 

The output state in Figure 1, which is in general dif¬ 
ferent from the input state due to the action of DUT, can 
then be written as 


P’ = <P). (2) 

where e is an operator representing the effect of the DUT 
on the input state. This can be further expanded as |2] 

e(p) = '^E^pE\, 
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( 3 ) 
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where Ei's comprise a set of at most four operators de¬ 
scribing the effect of DUT. Now these operational ele¬ 
ments can be expressed in terms of a fixed set of basis 
operators, Ek, fc = 1, 2, .A, i.e., we can write 


4 

Ei - ^ ^ eijYlEyyi, (4) 

m—1 


As a result, 


P — ^mpEj^Xmn} (5) 


where 


Xmn / , 


Since indices m and n each run from 1 through 4, Xmn 
is a 4 X 4 matrix, called the process matrix. This matrix 
is Hermitian = x)- Therefore, it has at most 4^ = 16 
independent parameters. Additionally, it is non-negative 
definite, i.e., its eigenvalues are zero or greater. 


Now, invoking the fact that for a trace-preserving pro¬ 
cess, Tr(p') = 1, one obtains from Eq. 5 

Y,XmnEiEm=I. ( 7 ) 

mn 

These are, in effect, four constraints on the elements, 
Xmn- These constraints then reduce the number of inde¬ 
pendent parameters of the x matrix from 16 to 12. In 
general, including nontrace-preserving processes Hi, 

P < I, (8) 

where 

P = 'y EiPi = y XmnEl^Em- (9) 

i mn 

Note that matrix P is Hermitian. 

In what follows, we choose Ei = di, where ui = J, 
<72 = <Jx, <73 = Uy, and cr4 = (7z, i.e., the identity matrix 
and the three Pauli spin matrices form the set of basis 
operators (this is a common choice in quantum comput¬ 
ing). It can be shown that for this fixed set of basis of 
operators Tr(x) = Ti"{P)/‘^i which then equals one for 
a trace-preserving process because in that case, P = I. 
Note that the trace of the x matrix may not always equal 
one as it depends on the choice of basis operators as well 
as the trace-preserving property of the quantum process. 

Now Eq. 8 implies that the eigenvalues of the P matrix 
(defined in Eq. 9) are each less than or equal to one. For 
the choice Ei = <Ji, we find 


Tr{x) + 2 V(/to(x 34) -b Re{xi 2 )Y + (/to(X 24) - Pe(xi3))^ + {Im{x 2 ?.) + ReixuW < 1, (10) 

rr(x) - 2 a/(/ m(x34) + ReixuW + (7to(X 24) - Pe(xi3))^ + (7m(x23) + i?e(xi4))^ < 1- (H) 


The left hand side of Eqs. 10 and 11 are the expres¬ 
sions for the eigenvalues of the P matrix in terms of the 
X matrix elements. Tr{x) = Xii + X 22 + X33 + X44- 
When the process is trace-preserving, the equality in 
Eqs. 10 and 11 holds, which then requires that all three 
terms under the radical sign be individually equal to 
zero because Tr(x) = 1- In other words, constraints 
i) /m(x34) = -Pe(xi2), h) 7 to(x24) = Reixis), and iii) 
/m(x23) = —7?e(xi4) must also hold in any numerical fit 
to the experimental data to yield a physical x matrix. 
These three constraints, along with the well-known con¬ 
straint, Tr(x) = I, comprise the four constraints of Eq. 
7 [9]- To our knowledge, this complete set of constraints 
in the above simple, explicit forms has not been cited or 
discussed in the past. Rather, sets of equations of the 
form, Eq. 7, seem to have been employed directly as 
constraints in numerical optimization procedures, using 
Lagrange multipliers, to obtain a fitted physical process 


matrix from experimental data (see, e.g., [lllin]). Knowl¬ 
edge of this full set of constraints, especially of the three 
additional explicit forms expressing relationships among 
the off-diagonal elements, provides additional insight into 
the structure of a process matrix. These explicit forms 
of the constraint equations can then serve as a useful 
tool in understanding and interpreting the experimental 
tomographic data. 

Along with being nonnegative Hermitian, the x matrix 
must, in general, satisfy the above two constraints. Eq. 
10, however, is a much stricter constraint than Eq. 11 
(assuming a positive sign for the term involving the rad¬ 
ical sign), which implies that the latter is redundant. In 
other words, Eq. 10 suffices. This is to be contrasted with 
the less stringent result Tr{x) < 1 [II], which is normally 
quoted in literature (see, e.g., [EldS]). Thus, in general, 
for a quantum process known to be nontrace-preserving 
like the polarizer (where Tr(x) = 1/2 ideally), or for a 
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process suspected to be not strictly trace-preserving like 
a quantum gate with leakage errors, or simply for a DUT 
whose behavior is not known a priori, the general con¬ 
straint, Eq. 10, must hold (or be explicitly invoked) in 
any fit to the experimental data. 

Recall that the process matrix for a trace-preserving 
process has 12 independent parameters, so 12 differ¬ 
ent measurements on the output quantum states are re¬ 
quired. For non-trace preserving processes, one must also 
include the effect of attenuation for each of the outputs 
giving rise to four additional parameters for a total of 
16. This is clear from considering the characterization of 
each of the four output states by quantum state tomogra¬ 
phy. For the purpose of discussion, let us consider opti¬ 
cal polarization qubits. In this case, a common choice 
is to prepare probe states of horizontal (|0)), vertical 
(|1)), diagonal ([|0) -I- |l)]/-\/2), and circular polarizations 
[|0) -I- i\l)\/^/2). QST of the single qubit requires four 
independent measurements (corresponding to the four 
classical Stokes parameters for each of the four output 
states). Thus a total of 16 measurements are needed to 
characterize the states output from the DUT. 

As is well known, in practice, due to experimental er¬ 
rors, the four initial states may not be prepared exactly 
as indicated. For example, the input horizontal may not 
exactly correspond to a pure |0) state, but in fact may be 
an admixture of |0) and |1) states. Such a mixture could 
either be coherent corresponding to an angular deviation 
from true horizontal, or it could be incoherent, corre¬ 
sponding to a slight depolarization error. Thus quan¬ 
tum state tomography should be performed on the input 
states as well as the output states. 

We now describe the procedure for the construction 
of the experimental x matrix using measured input and 
output states. It follows the procedure given in [2], but 
is modified by the fact that the input states are measured 
and the process may be nontrace-preserving. Let pj[j = 
1,2,3, and 4) denote the four experimentally prepared 
and measured input states. We now write 

Pj ( 12 ) 

i 

The Tji coefficients are then easily determined from 
Eq. 12. In this way, each measured input state is decom¬ 
posed into the basis we have chosen and will be included 
via TjiS in the calculation of the QPT below. 

The measured output states, expressed in terms of the 
fitted density matrices p", are first scaled to account for 
the possibility of a nontrace-preserving process (see the 
discussion pertaining to Eq. 1): 

j-{out) 


where / denotes the intensity of the qubit flux, and the 
superscripts {out) and {in) refer to the output and in¬ 
put states, respectively. In the quantum photon counting 


case, I is replaced with N^ where N denotes the number 
of photons counted per constant measurement interval. 
As in Eq. 12, we now write 

p'j ='^>^jkcrk, (14) 

k 

where the coefficients Xjk are determined from the above 
equation. Now inserting Ei = ai (our choice of basis 
operators) in Eq. 4, we have 

Pj = '^rnPjCriXmn- (15) 

mn 

Substituting Eqs. 12 and 14 into Eq. 15, we obtain m 

^ ^ ^jk^k — ^ ^ ^ ^ jiO'iO'n\mn ■ ( 1 ^) 

k mn i 

We now find we can write 

rji<Tm<JiCrn = PJT'crk, (17) 

i 

since Cfe’s form a complete set; /3^” are complex coeffi¬ 
cients. Substituting Eq. 17 in Eq. 16, and simplifying, 
we obtain 

(18) 

mn 

If we now regard as forming a 16x16 matrix, whose 
columns are indexed by mn and rows by jfc, then the 
above equation can be rewritten in a compact form as 

A = PZ (19) 

where A and x &re 16 x 1 column vectors whose rows are 
indexed by jk and mn, respectively. Inverting the above 
equation, we obtain 

X = kX, (20) 

where k = /3~^. Alternatively, we can write 

Xmn = ^ian^jk- (21) 

jk 

Having obtained x experimentally via the determination 
of the Tji parameters from Eq. 12, the A parameters from 
Eq. 14 and the /3 matrix from Eq. 17, it is critical to en¬ 
force the following required properties of the process ma¬ 
trix when fitting the experimentally determined process 
matrix (Eq. 21): 1) x is Hermitian, 2) it is nonnegative, 
and 3) its trace satisfies Eq. 10. 

The above requirements on the experimentally deter¬ 
mined process matrix can be enforced via conveniently 
available convex optimization software packages. After 
one has gone through the above “correction” procedure 
on the experimentally determined x matrix, one can ex¬ 
tract the operation elements, Ei, characteristic of the 
quantum device as shown in [5]. 
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We have revisited the theoretical aspects of single qubit 
quantum process tomography to determine the behavior 
of a quantum device. More specifically, we have reexam¬ 
ined the well-known constraints for the process matrix, 
and recast them into more insightful forms. In the case of 
a trace-preserving process, specific relationships among 
the various elements of the process matrix emerge that 
then shed light on its basic generic structure. Knowl¬ 
edge of these new constraint relationships permit an en¬ 
hanced understanding of the interpretation and analy¬ 


sis of the experimental data. Moreover, we point out 
a tighter bound on the trace of the process matrix for 
a nontrace-preserving process. Implementing the new, 
simplified form of constraints also serves as an alternative 
way to numerically fit the experimental data to obtain an 
optimal, physical process matrix. Additionally, for com¬ 
pleteness, we have described how to compute the pro¬ 
cess matrix with experimentally measured input states, 
which goes beyond the general assumption of perfect in¬ 
put state preparation within the framework of quantum 
process tomography (see, e.g., [5]). 


[1] I. L. Chuang and M. A. Nielsen, Journal of Modern Op¬ 
tics 44 , 2455 (1997) 

[2] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information: 10th Anniversary Edition, 
10th ed. (Cambridge University Press, New York, NY, 
USA, 2011). 

[3] K. Banaszek, G. Dariano, M. Paris, and M. Sacchi, Phys¬ 
ical Review A 61 , 010304 (1999). 

[4] D. F. V. James, P. G. Kwiat, W. J. Munro, and A. G. 
White, Phys. Rev. A 64 , 052312 (2001) 

[5] R. Bhandari, “On single qubit quantum state tomogra¬ 
phy,” (2014), arXiv: 1407.6668 [quant-ph] 

[6] N. Peters, J. Altepeter, E. Jeffrey, D. Branning, and 
P. Kwiat, Quantum Info. Comput. 3, 503 (2003) 

[7] J. Altepeter, E. Jeffrey, and P. Kwiat, in Advances In 
Atomic, Molecular, and Optical Physics, Vol. 52, edited 
by P. Berman and C. Lin (Academic Press, 2005) pp. 105 


- 159. 

[8] I. Bongioanni, L. Sansoni, F. Sciarrino, G. Vallone, and 
P. Mataloni, Phys. Rev. A 82 , 042307 (2010) 

[9] Note that these explicit forms of the four constraints 
could also have been derived directly by solving the four 
equations that Eq. 7 comprises. 

[10] J. L. O’Brien, G. J. Pryde, A. Gilchrist, D. F. V. James, 
N. K. Langford, T. C. Ralph, and A. G. White, Phys. 
Rev. Lett. 93, 080502 (2004) 

[11] This result is also obtained by simply adding Eqs. 10 and 

11 . 

[12] M. Mohseni and D. A. Lidar, Phys. Rev. Lett. 97, 170501 
(2006) 

[13] A. G. Kofman and A. N. Korotkov, Phys. Rev. A 80 , 
042103 (2009) 

[14] This expression and thus the definition of the /J coeffi¬ 
cients (see Eq. 17) are slightly different from the ones in 

m- 




